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MATHEMATICAL REPRESENTATION OF CURRENT DENSITY 
PROFILES FROM ION THRUSTERS 
byThaineW. Reynolds 
Lewis Research Center 

SUMMARY 

A concise and convenient way of expressing ion beam spreading characteristics is 

presented. The two-parameter equation, j(z, 0)/[j(z, 0)] = e"^ 1 " cos is shown to 

be capable of representing experimental beam current density profiles from several bom- 
bardment ion thrusters. This is a far -field expression. A calculation of the development 
of the profiles downstream of an extended source indicates that ion beam profiles should 
be taken at least eight thruster radii away to be representative of the far -field behavior 
of the beam. 

Parametric plots (involving various combinations of n and A.) of thrust factor, ion 
beam envelopes, and back -sputtering due to ion beam erosion illustrate the use of the ion 
beam equation in beam -interaction calculations. Experimental ion beams illustrated had 
thrust factors varying from 0.95 to near 1.0. 

INTRODUCTION 

A description of the ion beam density profile of electron bombardment thrusters is 
needed for both mission analysis and spacecraft design. The angular spread of the ion 
beam as it leaves the thruster affects placement of possible spacecraft components, such 
as solar cells and thermal control surfaces, to avoid the problems of the ion beam strik- 
ing sensitive surfaces and causing sputtering erosion. A loss in thrust is also involved 
when a significant fraction of the beam ions have trajectories away from the axis. 

This report presents a two -parameter equation which has been found useful in ex- 
pressing analytically the ion beam profile of a number of bombardment ion thrusters for 
which experimental profile data were available. The procedure for determining the 


values of the parameters for fitting the equation to experimental data is described. Fi- 
nally, several calculations involving interactions of the beam profile are presented to il- 
lustrate its utility. 


DENSITY PROFILE EQUATION 
Form of Equation 

Ion beam density profiles from electron bombardment thrusters have been described 
as having a relatively flat central core and an exponential dropoff with increasing distance 
from the beam centerline (ref. 1). An equation which exhibits these characteristics is 

j(z, e) = e -[x(i - cos e)f (j) 

j(z,o) 

where as shown in sketch (a), 



j(z , e) current density in radial plane at axial distance z and radial angle 6 (The cur- 
rent density is taken in the plane normal to the z -direction (radial plane) though 
the particle trajectories are at angles to this plane.) 

j(z,0) current density in radial plane at axial distance z on centerline 
6 half -angle from centerline to point of measurement 

n, X parameters to be determined 
(All symbols are defined in appendix A. ) 

This equation is deduced, then, simply from observation that experimentally deter- 
mined profiles can be reasonably well described with the appropriate choice of the param- 
eters n and Equation (1) is analogous to a similarly derived expression for density 
profiles of supersonic jets expanding into a vacuum (ref. 2). 
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Figure 1 shows the typical shapes of profiles obtained for n = 1 and 2 and for several 
values of X. A higher X for any given n indicates a more collimated beam. Similarly, 
a higher n for any given X indicates a more collimated beam. 

Equation (1) is a far -field expression. It describes the expanding flow as if the ion 
trajectories were straight lines which appear to originate from a point source. The cur- 
rent density thus would vary along any trajectory inversely as the square of the distance 
from the source. The flow is completely defined as a function of the angle Q. Equation 
(1) must necessarily become inapplicable close to the thruster. The distance from the 
thruster within which the equation is not useful is discussed in appendix B. 

The total current at any axial location z must be constant and equal to the emitted 
current. Thus, 




Substitution of equation (1) into equation (2) and rearranging yield 


Letting 


3 (z , 0) 

I 




e -[X(l - cos 


rtf r dr 


y = 


r 

z 


equation (3) becomes 



Designating the bracketed term in equation (5) as c(n, X), 


( 2 ) 


(3) 


(4) 


(5) 
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1 


c(n, X) = 




y dy 


( 6 ) 


equation (5) can be expressed 


j(z,0) 

j 



( 7 ) 


1 Finally, equations (7) and (1) combine to give the complete current density relation 
at any location, 


= c(n, X) ^2 e"^ 1 ~ cos (8) 

J z 2 . 

i n 

The part I of equation (8) gives the variation along the centerline; part n gives the ra- 
dial variation. Values of c(n, X) are plotted in figure 2. 


Comparison of Equation with Experimental Data 

The applicability of equation (1) to describing experimental ion beam profile data is 
shown in figures 3(a) to (g). In these figures the experimental profile data cited (refs. 1 
and 3 to 6) are shown by the data symbols. The solid lines are the plots of equation (1) 
for the particular values of n and X cited on each curve. The value of z/Rq, the ratio 
of the distance downstream that the experimental profile was taken to the thruster radius, 
is also shown. 

Two sets of experimental data for a given thruster have ion beam density profiles at 
two or more axial distances (figs. 3(d) and (e)). In figure 3(d), the data at axial distance 
ratios z/Rq of 14. 8 and 25. 2 fall on the same curve, showing that the assumption of lin- 
ear trajectories between the two distances was applicable to this case. In reference 1, 
from which the data of figure 3(e) were obtained, the thruster size is not given. The ac- 
tual linear distance at which the several profiles were taken is therefore shown in figure 
3(e), rather than the distance ratios z /Rq which are cited in the other plots. The pro- 
files at axial distances of 16. 5, 21. 5, 28.4, and 41.4 centimeters are closely approxi- 
mated by the single profile equation with n = 1. 19 and X = 85. The data of the closest 
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profile (z = 11. 5 cm) may be too close to the thruster to expect far -field behavior (see 
appendix B). The furthest downstream profile (at z = 53.4 cm) shows considerable asym- 
metry and is plotted separately in figure 3(e) along with the same equation plot as in the 
other part of this figure. The profiles of figure 3(e) also support the assumption of linear 
ion trajectories in the far field. Reasons for the asymmetry of the one profile were not 
discussed in reference 1. 

The other profiles of figure 3 were taken at distance ratios z/Rq varying from 4 
to 23.3. 

It can be seen that equation (1) is capable of representing the experimental, ion beam, 
current density profile data reasonably well. The greatest differences between experi- 
ment and equation values occur at the low -density values where the accuracy of determi- 
nation is least. 

Some of the information from figure 3 is summarized in table I. 


Determination of Parameter Values 

An experimental profile can be readily checked against equation (1) to determine the 
appropriate values of n and A. The value of A is determined, independently of n, from 
the angle e e at which j(z, S e )/[j(z, 0)] = ^(i.e., 0.368), 

A = (9) 

1 - cos e e 

Then, from any other experimental point, the value of n is determined. For instance, 
at the angle 0 Q j, at which j(z, e Q 1 )/(j(z, 0)) = 0. 1, 


ln[A(l - cos 9q j)] 

The radial profile will thus fit at the three values of angle 0 , 0, 9 e , and ^ . 

APPLICATION TO CERTAIN PROBLEMS 
Thrust Factor 

The ratio of the axial component of thrust to the thrust that would be obtained if all 
trajectories were directed along the axis is referred to herein as the thrust factor Frp. 
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Thus, 


F T = 


£ 


m(r)v cos 6 2nr dr 


/ OO 

m(r)v27ir dr 


(ID 


The velocity v is assumed to be constant for any trajectory and to be directed in the Q- 
direction. Equation (11) can be expressed, then, as 



( 12 ) 


Substituting equations (1) and (4) into equation (12) yields 


exp 


F T ~ 


XI - 


Vi^y 2j 




Vi + y 2 


dy 


exp < - 


'0 


|A 1 - 


Vi + y 2i 




(13) 


>Y dy 


Numerical integration of equation (13) resulted in the values of thrust factor shown plotted 
in figure 4. 


Ion Beam Envelope 


The fraction of the total ion beam enclosed within a given half -angle 0-^ is expressed 


as 


6 
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(14) 


J (g L ) 

I(total) 


l 


•r=z tan 0 L 


,-[A.(l - cos 6)f 2 „ ir 


f 


-[X(l - cos 8)f 2ur dr 


Substituting equations (4) and (6) into equation (14) yields 


I(total) 


= 2c (n, 



(15) 


Numerical integration of equation (15) for values of n = 1 and 2 and for several A -values 
yielded the curves shown in figures 5(a) and (b), for the fraction of the total beam con- 
tained within the angle 8-^. 


Return of Sputtered Target Flux to Thruster 

The testing of ion thrusters in vacuum tanks causes sputtering of the facility walls 
where the ion beam impinges. The return of this sputtered material to the thruster face 
and condensation thereon was shown to affect the operation of glass -coated accelerator 
grids (ref. 7). In order to get some estimate of the amount of back -sputtered material to 
expect, the following model was used (sketch (b)), which made use of the proposed beam 
profile equation (1): 
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The returning (back sputtered) flux rate p 2 depends on the local emission rate of 
sputtered material and the view factor dF^ ^ from the point of sputtering to the 
thruster center. (The variation across the thruster face will not be more than a few per- 
cent as long as the sputtered surface is about five or more thruster radii away. ) 


p 2 = 



i^dF 


1,2 


(16) 


Assuming cosine -law distribution of the sputtered material and cylindrical symmetry, 
the back flux from the target and the wall can be written separately and summed. (View 
factors defined in appendix A are derivable from relations in reference 8. ) 



The sputtered emission rate depends on the incident ion rate ^ and the sput- 
tering ratio S, 


= MjS ( 18 ) 

In turn, S depends on incident ion energy, target material, and angle of incidence. 
The angular variation was assumed to follow the relation (ref. 9) 

S = S(0) sec a (19) 

where a is the angle between the surface normal and the incident ion direction. 

For the variation of with the position variables z and r, use was made of 
equation (8). 

It can be seen in sketch (a) that for the circular target a = 0, and that for the cylin- 
drical tank wall a = ( n/2 ) - 0. Also the arrival rate on the cylindrical wall is related to 
the arrival rate on the target at the same location by 

k ,R)l = [mi(z,R) 1 tan 9 (20) 

J Tank wall L 1 JTarget 
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By converting to the dimensionless variable x, where 


x = 


z 

R 


( 21 ) 


and using relations (4), (8), and (18) to (21), equation (17) can be transformed into 


l±2 

~v 








+ 2S(0) 

V. 



( 22 ) 


Term I of equation (22) is the back flux from the circular target. Term II is the back 
flux from the cylindrical wall. The two parts of equation (22) are plotted separately in 
figures 6 and 7 for values of n = 1 and 2 and for several values of A. 

In figures 6(a) and (b), for example, the ordinate parameter 


v 2 /M 2 1 

v \R 0 / 2S(0) 

is plotted. This parameter accounts for the flux ratio \i^/ v > f rom a circular target of 
radius r to the thruster exit center. 

The return flux from a target is obtained, then, as 



target 


^2 


J T 1 

r 2 2[S(0) ] 
0 


target 


^2 2S(0) 

4 


(23) 
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It can be seen, for instance, that for n = 1 and A = 50, no significant fraction of the 
total return flux from a circular target is obtained beyond about a 20° half -angle. For 
n = 1 and A = 15, most of the return flux is confined to the target that is within a 6 -value 
of about 45°. For n = 2 and these same A-values, the cutoff angles are nearer 15° and 
30°, respectively. 

Figures 7(a) and (b) show the return flux from the tank walls to the thruster exit cen- 
ter. Here the ordinate gives the return flux from x = 0 out to any axial distance x = z/R 
(The reference radius R in these plots is the tank radius, not the thruster radius. ) The 
return flux from the cylindrical wall is then 



wall *- 


l M .2.ir _ 

7 r2 2S «»| 


2S(0) 


wall 


R 


(24) 


In this case it can be seen that for n = 1 and A = 50 there is essentially no return 
flux from the wall between x = 0 and x = 1.5. The major portion of the return flux from 
the wall comes from the wall area between x = 1. 5 and x = 10. For n = 1 and A = 15, 
return flux from the wall close to x = 0 is noted, but essentially no wall -return flux 
comes from the wall beyond about x = 8. For n = 2, the same A-values represent more 
collimated beams; hence, there is less wall -return flux from similar facility dimensions. 


CONCLUDING REMARKS 

A two -parameter equation is shown to be capable of representing experimental ion 
beam density profiles of bombardment ion thruster exhausts. The equation is a far -field 
relation (eq. (1)) and therefore will not be applicable too close to the thruster exit. A 
calculation of the development of the profiles downstream of an extended source indicates 
that ion beam density profiles should be taken at least eight thruster radii away to be rep- 
resentative of the far -field behavior of the ion beam. The equation was shown to repre- 
sent profiles taken as close as four thruster radii downstream. However, these may not 
be representative of the far -field patterns. 

The particular advantage of being able to express the ion beam profile in a concise 
analytical relation has been illustrated by parametric plots of thrust factor and ion beam 
envelope curves. Experimental beam profiles illustrated had thrust factors varying from 
about 0. 95 to nearly 1. 0. 
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Use of the beam profile equation to calculate back -sputtering from vacuum chamber 
walls during ion thruster testing was illustrated. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 18, 1971, 

129-02. 
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APPENDIX A 


c(n, X) 



I 

i(e) 

I(total) 
}(z, 9) 


j(z,r) 


i(z,o) 

j(0,r Q ) 

3 ( 0 , 0 ) 

f 

m(r) 

n 

R 

R 0 

r 

r 0 

s 

S(0) 

v 

x 

y 


SYMBOLS 

parameter defined in eq. (6) 

o 

differential area of thruster, cm 

differential view factor, fraction of emitted flux from surface 1 that inter- 
cepts surface 2 

current, A 

current enclosed within cone of half -angle 0, A 
total current, A 

current density in radial plane at axial distance z and radial angle 0, 
A/cm 2 

current density in radial plane at axial distance z and radial distance r, 
A/cm 2 

p 

current density in radial plane at axial distance z on centerline, A/cm 
current density in the radial plane at thruster exit (z = 0) and at radial 

9 

distance Tq, A/cm 

current density in radial plane at thruster exit (z = 0) at centerline, A/cm 2 

O 

average current density of thruster (eq. (B3)), A/cm 

p 

mass flux in radial plane at radius r, kg/(cm )(sec) 
parameter of eq. (1) 

radius of facility vacuum tank (sketch (a)), cm 
radius of thruster, cm 
radius, cm 

radius in thruster exit plane, cm 
sputtering yield, atoms/incident ion 

sputtering yield at normal incidence, atoms/incident ion 
velocity 

nondimensional axial distance, z/R 
nondimensional radial distance, r/z (= tan 0) 
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Of 

J8 

e 


d e 

° 0 . 1 
0 L 


X 

Mi 

M 2 

*1 


axial distance variable, cm 

distance of target from thruster exit plane, cm 

angle between surface normal and incident ion direction 

angle in thruster exit plane (sketch (c)) 

angle from centerline of ion beam 

angle at which j(z, 0 e )/j(z, 0) equals 0. 368 

angle at which j(z, e Q 1 /j(z, 0) equals 0. 1 

cone half -angle enclosing given fraction of total ion beam 
(eqs. (14) and (15)) 

parameter of eq. (1) 

particle arrival rate at surface in target plane, atoms/(cm )(sec) 

p 

particle arrival rate at center of thruster exit plane, atoms/(cm )(sec) 

p 

average emitted ion rate from thruster, ions/(cm )(sec) 

p 

particle emission rate from target surface, atoms/(cm )(sec) 
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APPENDIX B 


DEVELOPMENT OF FAR -FIELD PROFILE FROM EXTENDED SOURCE 

In order to get some indication of the closest approach to an extended source (the 
thruster) that a far -field expression such as equation (1) could be useful, a calculation of 
the profiles as a function of distance from such a source was carried out. 

If it is assumed that nothing influences the ion trajectories once they are beyond the 
neutralization plane of the thruster, the development of a particular far -field profile, ex- 
pressible in the form of equation (1), would only come about because the local emission 
sources (for instance, the individual grid holes) exhibit this same type of profile. The 
amount of emission current from the individual sources may vary. In this particular ex- 
ample a distribution of current density across the radius of the thruster exit (z = 0) is as- 
sumed to be 


(Bl) 

1(0,0) \2 RqJ 

That is, individual emitting holes are not assumed, but rather a continuous emission 
source over the thruster face is assumed. The magnitude of the emission current varies 
as equation (Bl) and the distribution locally varies as equation (1). 

For the profile of equation (Bl) the current density at the centerline j(0, 0) is related 
to the average current density j by 


3(0,0) = 2. 16 j 


where 


» R o 

A schematic showing the geometric relations is shown in sketch (c). 


(B2) 


(B3) 
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The flux to point z, r from the differential area of the thruster is 

= C (n, A) ^2 e-W 1 - cos (B4) 

W’V ttz 2 


By using 


COS O' = 


z 



+ r 


- 2r Q r cos /3 + z 


2 


and 


j (0, r Q ) = 2. 16 j cos 



(B5) 


(B6) 


equation (B4) can be written 


j(z,r) _ 2. 16 c(n, A) 


] 


7TZ 



in 


+ r 2 - 2^r cos /3 + z 2 ^ 


r 0 dr 0 d /3 


(B6) 


Equation (B7) was integrated numerically. Results for n = 1 and X = 20 are shown 
plotted in figures 8 and 9. Figure 8 shows the variation of j(z,r)/j with distance for sev- 
eral e -values. Values of j(z,r)/j approach closely those given by the far -field expres- 
sion within about six to eight thruster radii downstream. The same calculations are pre- 
sented in a different manner in figure 9. This figure shows constant current density lines. 
The near -field data (dashed lines) blend into the far -field values (solid lines) by about six 
to eight thruster radii downstream. 
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TABLE I. - COMPARISON OF ION BEAM PROFILES WITH EQUATION (1) 



Thruster identification 


Axial distance 

Constants 

Thrust 

Diameter, 

Type 

Figure 

Refer - 

of profile, 

z/R Q 

n 

A 

factor , 
Ft 

cm 



ence 





30 

Glass grid 

3(a) 

3 

10 

1 

24.3 

0.954 

5 

Glass grid 

3(b) 

4 

23. 3 

.92 

40. 2 

.966 

5 

Glass grid 

3(c) 

(a) 

4 

1.28 

18. 2 

.951 

5 

Two grid 

3(d) 

5 

14. 8 

.9 

63. 3 

.980 





25. 2 

.9 

63.3 

.980 

(b) 

(b) 

3(e) 

1 

— 

1. 19 

85 

.99 

15 

Two -grid 

3(f) 

6 

4 

1. 23 

34.6 

c . 976 






1.4 

17.0 

d . 955 

30 

Glass grid 

3(g) 

7 

4. 8 

. 8 

42.9 

.960 


a Unpublished data by J. H. Molitor and H. J. King of Hughes Research Laboratory, Malibu, Calif. 


^Not given in reference. 
c Case HI of ref. 6. 
d Case I of ref. 6. 
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Current density ratio, j(z, 0)/j(z, 0) 




Current density ratio, j(z f 0)/j(z, 0) 



(a) Data of reference 3. Glass-grid, 30- (b) Data from reference 4. Glass-grid, 5- 

centimeter-diameter thruster; n = 1; centimeter-diameter thruster; n = 0. 92; 

A = 24.3; z/Rq =10. A - 40. 2; z/R 0 = 23.3; experimental profile 

taken 58. 4 centimeters downstream. (Data 
from fig. 10 of ref. 4). 



Angle from centerline of ion beam, 0, deg 

(c) Unpublished data by J. H. Molitor and (d) Data from reference 5. Glass-grid, 5- 

H. J. King of Hughes Research Laboratory, ' centimeter-diameter thruster; n = 0.9; 

Malibu, Calif. Glass-grid, 5-centimeter- A =63. 3. 

diameter thruster; n - 1.28; A - 18. 2. 

Figure 3. - Comparison of profile equation (1) with data. Solid line is equation for n, A values cited. 
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Angle from centerline of ion beam, 0, deg 


(e-1) For various axial distance z. (e-2) For axial distance 

of 53.4 centimeters. 


(e) Data of reference 1. Thruster type and diameter not specified; 
n a 1.19; \ c 85. 



Angle from centerline of ion beam, 0, deg 


(f) Data of reference 6. Two-grid, 15- (g) Unpublished data by R. Bechtel 

centimeter-diameter thruster; z/Rn = 4.0. of Lewis. Glass-grid, 30- 

centi meter -diameter thruster; 
n » 0. 8; X * 42. 9; z/R 0 = 4.8; 
thrust factor, Fy = 0.960. 

Figure 3. - Concluded. 
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Dimensionless radial distance, r/R n Current density ratio, j(z,0)/} 



Figure 8. - Variation of current density with distance at several angles. Straight lines are far-field values; 
data points are calculated values from extended source relation (eq. (B7)); n * 1; A= 20. 



0 2 4 6 8 10 12 14 16 

Dimensionless axial distance, z/Rg 


Figure 9. - Constant current density contours. Solid lines from equation (1); dashed lines from equation 
(B7); n= 1; A= 20. 
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